Entropy scaling and thermalization in hadron-hadron collisions at LHC 
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The recent LHC data have once more brought the issue of the formation of thermalized state of 
matter in hadron-hadron collisions into the forefront. In this letter, we have shown that the scaling 
of the information entropy of the chaotically produced particles is valid up to recently available 
data from p+p collisions at -^s = 2.36 TeV obtained by ALICE experiment. We predict that at the 
highest energies projected for LHC, almost all the particles will be produced chaotically, indicating 
that a collective behavior should be observed in hadronic collisions, as indicated by one of present 
authors quite some time ago. 
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The properties of strongly interacting matter at ex- 
treme conditions have been an intense area of research 
for quite some time now. The results from relativistic 
heavy-ion collision experiments [J have enriched our 
knowledge of matter at high temperatures and small den- 
sities. One of the main goals of these experiments is to 
test our understanding of Quantum Chromodynamics, 
the fundamental theory and look for the formation of 
deconfined quark-gluon matter. 

The non-perturbative nature of QCD makes it diffi- 
cult to understand all aspects of particle production in 
strongly interacting systems. The huge number of parti- 
cles produced in heavy-ion collision experiments makes it 
even more difficult to decipher the proper sources of these 
particles. Under such circumstances, one can look at the 
simpler systems, such as hadron-hadron collision to gain 
better understanding of particle production mechanisms 
and their possible sources. 

Though the research on the multi-particle productions 
in different systems, such as high energy hadron-hadron, 
hadron-nucleus and nucleus-nucleus collisions are being 
pursued by various groups for a long time (see e.g. [3| 
and references therein) , a proper theory of multiparticle 
production remains elusive till date. More so, as non- 
perturbative nature of the interactions involved makes it 
difficult to calculate the soft processes. One of the pos- 
sible routes to the understanding of multiplicity is to in- 
vestigate the equation of state of the matter produced in 
these collisions [1,[5|. This, along with a hydrodynamic 
picture to take care of the overall energy-momemtum 
conservation during the space-time evolution of the ex- 
cited matter, can give some idea of the dependence of 
multiplicity on energy. 

Multiplicity distribution of hadrons produced in high 
energy particle collisions has long been known to devi- 
ate from a Poisson distribution and has been regarded as 
a potentially useful source of information about the un- 



derlying production processes. Multiplicity distributions 
have been predicted using various models of the low -pt 
hadron production processes in hadron-hadron collisions 

The Feynman scaling hypothesis led to the idea of 
KNO scaling ll|. However it may be noted that the 



inelastic hadronic interactions do not strictly follow the 
KNO scaling [13] ■ Inelastic cross section is usually 
made up of single diffraction dissociation and the non- 
single diffractive (NSD) part, the latter being the main 
source of particles produced. The collision energy de- 
pendence of the NSD multiplicity was nicely explained 
by this scaling up to the ISR energies. But this scaling 
law broke down 14 , once the collision energy reached 
the higher SPS energies. 

Analysis of the multiplicity data using statistical mo- 
ments \l4 . [lit ^^-d scaling laws provide new ideas to 
interpret the results. A ne w q uantity, the entropy of 
multiplicity distributions 16j, llTI ] , was proposed to revive 
the scaling law in high energy hadronic collisions at the 
collider energies. An important observation was the de- 
pendence of the multiplicity distribution on energy and 
rapidity bins based on the assumption of the existence of 
two kinds of sources 13, HI • Firstly, the chaotic source, 
which is concentrated at small rapidities and has the 
characteristics of thermally equilibrated system and sec- 
ondly, the coherent source, that contributes to the whole 
rapidity region. It was also shown that the new scaling, 
based on entropy, holds good for the entire range of en- 
ergy starting from ISR {^/s = 19 GeV) up to the highest 
available collision energy at SPS (-^/s = 900 GeV) only if 
the entropy is calculated from the chaotically produced 
particles [19|. With increasing energy, the weight of the 
chaotic or thermally equilibrated source is expected to 
increase. This may also be indicative of the production 
of QGP, one of the main goal of all these investigations. 
Recently, deconfinement in hadronic collisions at TeV en- 
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ergies has been explored using 1-d hydrodynamic model 
by E735 collaboration 

In this letter, we will extend the results from this scal- 
ing law [l9| for the NSD multiplicity distributions in p- 
p collisions at the highest available collision energies at 
LHC and try to extrapolate to the highest energy pro- 
jected at the LHC. We shall argue that even in hadron- 
hadron collisions, one can expect a thermalised state of 
matter at sufficiently high collision energies, as was al- 
ready alluded to in [1, Q . 

The scaling variable, information entropy has been cal- 
culated within the context of the two component model 
[itI from multiplicity distribution. According to this 
model the emission of particles occurs from a convolu- 
tion of a chaotic source with Planck-Polya distribution 
and a coherent source with Poisson distribution. The en- 
tropy for symmetric rapidity intervals jryj < rjc , for some 
arbitrary rjc, and at center of mass energy -^s is given by. 



- (?ich(??c,ys)) ln(nch(?7c, Vs)) (1) 

where rich is the chaotic fraction of the total multiplicity. 
(Similarly, Hco and ntot denote the coherent and the total 
multiplicity, respectively, so that ritot = n-ch + tIco-) nch 
may be obtained as a product of the total multiplicity 
with the chaotic fraction P as. 



{rich) = P{ritot) 
where the chaotic fraction is given by, 

P=[C2-{l + l/{ntot))Y'^ 



(2) 



(3) 



and C2 — {ritot^) I {ritot}'^ is the second moment of mul- 
tiplicity distribution. A scaling law then emerges if we 
plot the scaled entropy S/rimax as a function of the scaled 
variable ^ = r]c/ rjmax where 



= ln[(v^- 2m^)/mn] 



(4) 



In Fig. [T]we have plotted S/r]max as a function off for 
different collision energies starting from SPS IJ, |l5| to 
the highest available energy at LHC [213 . As can be seen, 
all the experimental data follow a single curve within 
the estimated statistical errors. A single scaling law is 
thus able to describe multi-particle productions up to the 
highest energy available for hadron-hadron collisions. 

Expecting that the scaling law holds at even higher en- 
ergies, an important consequence is that we can predict 
both the average multiplicity as well as the nature of the 
particles produced at higher energies. For this one has 
to extrapolate the available data to the desired -y/s for 
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FIG. 1: Entropy scaling from two component model. 

both the average total multiplicity (ritot) and the aver- 
age chaotic multiplicity (rich)- For the extrapolation in 
yG, we choose those three values of f for which experi- 
mental data are available for the highest ^/s—2.36 TeV at 
ALICE. Experimental data are not available at exactly 
these values of ^ at lower values of ^/s. They are first ob- 
tained through interpolation from the experimental data 
for (ritot) and the model data for (rich)- 
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FIG. 2: (rich) as a function of ^. The squares and circles 
correspond to data and interpolated points respectively and 
the bands show the limits of statistical errors. The dotted 
line is from the fit described in the text. 

In Fig. [2] we have plotted the (ric/i) as a function of ^ 
for different collision energies. A good fit is obtained for 
different collision energies with the logistic function. 



b — a 



{nch) 



1 _|_ e-(?-c)/d 



(5) 



where a, b, c and d are parameters to be fitted. With 
this fit, we could obtain the value of (Uch) at any value 
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of ^ for a given ^/s. Thus we could obtain the {rich) at 
the desired values of ^. 
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FIG. 3: {rich) as a function of collision energy. The squares 
and circles correspond to data and predicted points respec- 
tively and the bands show the limits of statistical errors. The 
dotted line is from the fit described in text. 

We can now plot the mean multiplicity of chaotically 
produced particles {rich), as a function of collision energy 
for different values of ^ (Fig. [3]) obtained from eqn. [5] 
and ALICE data at •\/i=2.36 TeV. Once the new scahng 
law for entropy has been established it follows that {rich) 
behaves as a power law with ^/s for large values of {rich)- 
Indeed a good fit is obtained with the formula, 
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FIG. 4: (ntot) as a function of collision energy. The squares 
and circles correspond to data and predicted points respec- 
tively and the bands show the limits of statistical errors. The 
dotted line is from the fit described in text. 





Parameters 


Correlations 


0.051362 


b = 0.24(3) 
c = 0.37(2) 


(6,c) = -0.991 


0.102724 


b = 0.53(4) 
c = 0.36(1) 


(&,c) = -0.991 


0.133542 


6 = 0.8(1) 
c = 0.35(2) 


(&,c) = -0.995 



TABLE II: Fit parameters of ntot with ^/s for different ^ 



{rich) = bV~s ' (6) 

where h and c are parameters to be fitted. The fit param- 
eters are given in the Table H] Using this fitting function 
one can now predict the value of {rich) at higher collision 
energies that are expected to be obtained at LHC in near 
future. 





Parameters 


Correlations 


0.051362 


b = 0.10(1) 
c = 0.47(2) 


(b,c) = -0.991 


0.102724 


b = 0.24(2) 
c = 0.45(1) 


(6,c) = -0.990 


0.133542 


b = 0.35(5) 
c = 0.43(2) 


(b,c) = -0.994 



TABLE I: Fit parameters of {rich) with ^/s for different ^ 

A similar exercise for the average total multiplicity 
{ntot) (Figure. gives the collision energy dependence 
of the average total multiplicity for different ^. We again 
find the best fit to have a similar functional form as in 
eqn. [HI The fit parameters are shown in Table |lTl 



The predicted values of both {rich) and {ntot) are 
marked on the respective figures Fig. [3] and Fig. |4l and 
listed in Table IIIII The table also shows the contribution 
of the multiplicity of the coherent source (rico)- At these 
high energies the {nco) is consistent with zero. The na- 
ture of {nco) is shown in Fig. [5j We see that the {nco), 
after an increase at the collision energies up to 900 GeV, 
indeed shows a dip towards higher energies. The four 
points are not enough to obtain a suitable functional form 
to fit these data. The line is just a guide to the eye. 

This analysis suggests that at the higher energies the 
chaotic part of multiplicity completely dominates over 
the coherent fraction. Data at higher collision energies 
that would be available shortly would test these prdic- 
tions. 

To summarise, we have established the scaling law 
for the information entropy obtained in a two source 
model, as a function of ^/s upto the maximum possible 
range available at present. The multiplicity of chaoti- 
cally produced particles has been evaluated from a two- 
source model in hadronic collisions starting from lowest 
SPS energy to the highest available LHC energy. Our 
study leads to the prediction that at the highest energies 
projected for LHC almost all the particles are produced 
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FIG. 5: (rico) as a function of collision energy. Dotted line is 
to guide the eye. 



(GeV) 




{ritot) 


{rich) 


{rico) 




0.051362 


6.7(2) 


6.6(2) 


0.1(3) 


7000 


0.102724 


13.4(4) 


12.5(4) 


0.9(6) 




0.133542 


17(1) 


16(1) 


1(1) 




0.051362 


8.7(4) 


9.1(4) 


-0.4(6) 


14000 


0.102724 


17.3(6) 


17.0(7) 


0.3(9) 




0.133542 


22(2) 


21(2) 


1(1) 



TABLE III: Different multiplicities with varying ^ for different 



chaotically i.e. the chaoticity fraction P, reaches almost 
unity, indicating the formation of a thermalized source 
and the possibility of observation of collective behaviour 
in hadronic collisions. 
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